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Abstract 

We review the impact of nuclear forces on matter at neutron-rich ex¬ 
tremes. Recent results have shown that neutron-rich nuclei become 
increasingly sensitive to three-nucleon forces, which are at the forefront 
of theoretical developments based on effective field theories of quantum 
chromodynamics. This includes the formation of shell structure, the 
spectroscopy of exotic nuclei, and the location of the neutron dripline. 
Nuclear forces also constrain the properties of neutron-rich matter, in¬ 
cluding the neutron skin, the symmetry energy, and the structure of 
neutron stars. We first review our understanding of three-nucleon forces 
and show how chiral effective field theory makes unique predictions for 
many-body forces. Then, we survey results with three-nucleon forces 
in neutron-rich oxygen and calcium isotopes and neutron-rich matter, 
which have been explored with a range of many-body methods. Three- 
nucleon forces therefore provide an exciting link between theoretical, 
experimental and observational nuclear physics frontiers. 
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1. INTRODUCTION 

Recent studies have shown that three-nucleon (3N) forces play a key role for understanding 
and predicting neutron-rich nuclei and for the formation and evolution of shell structure. In 
addition, the location of the neutron dripline, where nuclei cease to be bound, is sensitive to 
small interaction contributions and therefore to 3N forces. At the same time, 3N forces are 
the dominant uncertainty in constraining the properties of neutron-rich matter at nuclear 
densities, which is important for the structure of neutron stars. This leads to an exciting 
connection of 3N forces with the exploration of extreme neutron-rich nuclei at rare isotope 
beam facilities and with forefront observations in astrophysics. 

In the following, we discuss our understanding of nuclear forces based on chiral effective 
field theory (EFT) and show how this framework makes unique predictions for many-body 
forces. In particular, the properties of all nuclei can be predicted up to high order (next-to- 
next-to-next-to-leading order, N^LO), with just two new low-energy couplings in many-body 
forces. To this order all other three- and four-nucleon (4N) interactions depend on couplings 
that are determined from interactions with pions or enter in two-nucleon (NN) interactions. 
We discuss theoretical developments of the subleading 3N forces and advances towards 
including them in few- and many-body calculations. In Sectionsandwe survey results 
with 3N forces in neutron-rich oxygen and calcium isotopes, which have been explored with 
a range of many-body methods. These present key regions in the study of the neutron 
dripline, of shell structure, and in the spectroscopy of exotic nuclei. In Section]^ we discuss 
the impact of nuclear forces on the properties of neutron-rich matter, including the neutron 
skin, the symmetry energy, and the properties of neutron stars. Finally, we conclude and 
give an outlook with open problems and opportunities in Section]^ 
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Figure 1 

Chiral EFT for nuclear forces, where the different contributions at successive orders are shown 
diagrammatically [figure adapted from dEH. Nucleons and pions are represented by solid and 
dashed lines, respectively. Many-body forces are highlighted in orange including the year they 
were derived [3N forces at N^LO JsHlJ and N^LO JUIbll; 4N forces at N^LO (O]. All N^LO 3N 
and 4N forces are predicted parameter-free. 


1.1. Chiral EFT for nuclear forces 

Chiral EFT provides a systematic basis for strong interactions at momentum scales of 
the order of the pion mass Q ~ based on the symmetries of quantum chromodynamics 
(QCD) ([T][5|. In chiral EFT, nucleons interact via pion exchanges and shorter-range contact 
interactions. The resulting nuclear forces and consistent electroweak operators are organized 
in a systematic expansion in powers of Q/Ab, where Ab ~ 500 MeV denotes the breakdown 
scale, leading to a typical expansion parameter Q/Ab ~ 1/3. The EFT enables controlled 
calculations with theoretical error estimates, which is especially important for exotic nuclei 
and neutron-rich matter under extreme conditions in astrophysics. Moreover, chiral EFT 
connects nuclear forces to the underlying theory through lattice QCD (mini). 

Generally, nuclear forces are not observable and depend on a resolution scale A, so 
that the nuclear Hamiltonian is given by H{A) = T(A) -|- Vnn(A) -|- Vi 3 N(A) -|- V 4 n(A) ... 
As shown in Figure at a given order, nuclear forces include contributions from one- or 
multi-pion exchanges that govern the long- and intermediate-range parts and from short- 
range contact interactions. For each A, the scale-dependent short- range couplings are fit 
to low-energy data and thus capture all short-range effects relevant at low energies. While 
3N forces are not observable, there are natural sizes to many-body-force contributions that 
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Figure 2 

Different topologies that contribute to 3N forces up to N®LO. The shaded vertices denote the 
amplitudes of the corresponding pion/nucleon interactions. The individual diagrams are: (a) 27r 
exchange, (b) Irr-contact, (c) 3N contact, (d) 27r-l7r exchange, (e) ring contributions, 

(f) 27r-contact, and (g) relativistic corrections. Figure taken from Reference (I16II . 


are made manifest in the EFT power counting and which explain the phenomenological 
hierarchy of many-body forces, 14 jn(A) > V 3 n(A) > V 4 n(A) (1T][2|. The effects discussed in 
this review are dominated by the long-range parts of 3N forces and are therefore expected 
to be robust. Even though it is tempting to neglect contributions from 3N forces in cases 
when calculations based on only NN forces provide a good description of experimental data 
[see, e.g., Reference m], EFT power counting dictates the inclusion of all many-body 
forces up to a given order. In fact, explicit calculations show that 3N forces always provide 
important contributions in nuclei and matter (see Sections]^ [^and|^. The scale dependence 
can also be exploited by using the renormalization group (RG) to systematically change the 
resolution scale, while preserving low-energy observables. This can be advantageous for 
calculations of nuclei and nucleonic matter, because the evolution to lower scales facilitates 
the solution of the nuclear many-body problem due to a decoupling of low and high momenta 
in the Hamiltonian (Ell EH). 

1.2. Many-body forces 

Chiral EFT opens up a systematic path to investigate many-body forces and their impact 
on few- and many-body systems m- An important feature of chiral EFT is the consistency 
of NN and 3N (and higher-body) interactions. This determines the long-range two-pion- 
exchange parts of 3N forces at next-to-next-to-leading order (N^LO), with pion-nucleon 
couplings ci,C 3 ,C 4 , leaving only two low-energy couplings cd and ce that encode pion 
interactions with short-range NN pairs and short-range three-body physics, respectively dSl 
|4])- To fit CD and ce, different uncorrelated observables are used, e.g., the binding energy 
and half life of (1141) . or the binding energy of and the charge radius of '*He EH- 

At the next order, all many-body interactions are predicted parameter-free with many 
new structures, as shown in Figure These also depend on the leading NN contact inter¬ 
actions Cs,Ct (IH [6]) [see the 27r-contact contributions (f) and the relativistic corrections 
(g)]. Interestingly, for systems of only neutrons, the N^LO cd and ce parts do not con¬ 
tribute because of the Pauli principle and the coupling of pions to spin EH- Therefore, 
chiral EFT predicts all three-neutron and four-neutron forces to N^LO. 

The leading N^LO 3N forces improve few-body scattering, but interesting open problems 
remain EH). This makes the application of 3N and 4N forces at the next order (N^LO) very 
exciting. The derivation of N^LO 3N forces has only been completed recently (O 16]) , but 
no calculation exists for nuclei beyond A = 3, where the current state-of-the-art are N^LO 
NN plus N^LO 3N calculations, or consistent calculations at N^LO (see Sections and [^. 
The practical calculation of all topologies of Figure in a form suitable for few- and many- 
body calculations is a nontrivial task. Due to the large amounts of required computational 
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resources, 3N matrix elements so far have only been available in small basis spaces m- 
Recently, a novel and numerically much more efficient method has been developed for 
decomposing 3N forces in a plane-wave partial-wave basis GSJ. The new framework makes 
explicit use of the fact that all unregularized contributions to chiral 3N forces are either 
local, i.e., they depend only on momentum transfers, or they contain only polynomial non¬ 
local terms. These new developments allow to calculate matrix elements of all N^LO 3N 
contributions for large basis spaces required for ab initio studies of nuclei and nucleonic 
matter. An important open problem is to fit all NN and 3N low-energy couplings at N®LO 
consistently. To this end, it may be beneficial to fit several different few-body observables 
simultaneously within theoretical uncertainties, or to include also information beyond few- 
nucleon systems in the fits. The exploration of these strategies is currently in progress. 


2. NEUTROISI-RICH OXYGEN ISOTOPES 


With a closed Z = 8 proton shell, the oxygen isotopes provide an exciting laboratory to 
study nuclear forces with a range of many-body methods. The oxygen chain also exhibits 
remarkable aspects of exotic nuclei. In oxygen, the location of the neutron dripline at 
neutron number = 16 is anomalously close to the stable nuclei [at the same N as for 
carbon {Z = 6) and nitrogen {Z = 7)], see Section [2.Ij Adding a single proton to oxygen 
is able to bind six more neutrons to As discussed in Section 2.2 oxygen features two 
more shell closures in and ^'*0, where the latter doubly magic nucleus lies just at the 
neutron dripline. In recent years, experimental advances have even allowed explorations of 
oxygen isotopes beyond the dripline to ^®0, see Section 


2.3 


In the following, we discuss results for the neutron-rich oxygen isotopes and neighboring 
nuclei based on chiral EFT interactions. In particular, the softening of nuclear forces by RG 
transformations dniiTii) enabled many new calculations of nuclei with a range of many-body 
methods. The results with 3N forces surveyed in this review have been mainly performed 
using two different strategies for the Hamiltonian with different RG approaches. The first 
uses the RG to evolve NN interactions to low-momentum interactions Vjow k (Hill . At the NN 
level, the results discussed here start from the 500 MeV N®LO potential of Reference (PUjl . 
Ghiral 3N forces provide a general low-momentum basis, so that the low-momentum NN 
interactions have been combined with N^LO 3N forces fit to the binding energy and 
the "^He charge radius CSI. We will refer to these Hamiltonians as low-momentum NN-1-3N 
forces. Note that similar results are obtained if the RG evolution at the NN level is replaced 
by a similarity RG (SRG) evolution (I15II . The second strategy for the Hamiltonian uses a 
SRG transformation applied consistently to chiral NN and 3N interactions dUlllll). The 
results discussed in this review start from the same NN potential and include the same 
N^LO 3N operators, but with local regulators (f23ll. fit to the binding energy and half life 
of ®H (|14|l . We will refer to these Hamiltonians as SRG-evolved NN-|-3N-full Hamiltonian, 
or NN-|-3N-ind Hamiltonian if the SRG evolution starts only from the NN potential, but 
induced (ind) 3N interactions are kept. The latter corresponds to the results for the “bare” 
500 MeV N®LO potential, if induced higher-body interactions can be neglected. 

Next, we briefly discuss the different ab initio many-body methods used for medium- 
mass nuclei. These can be grouped into calculations based on valence-space Hamiltonians 
and calculations obtained directly in large many-body spaces. The latter treat all nucleons 
as active in a large basis space and rely on different, controlled approximations to solving the 
many-body problem of A nucleons. For calculations based on valence-space Hamiltonians, 
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the number of active degrees of freedom is reduced by treating the nucleus as a many-body 
system comprised of a closed-shell core with the additional valence nucleons occupying 
a truncated valence space. The valence-space Hamiltonian, which is diagonalized exactly, 
includes configurations from the large basis space via different many-body methods. 

Great advances to access large basis spaces have been made with coupled-cluster (CC) 
theory The CC method starts from a closed-shell reference state and includes corre¬ 
lations through a similarity transformation H = He^, where H is the normal-ordered 
Hamiltonian. In state-of-the-art calculations the cluster operator T = Ti + T 2 + ■ ■ ■ + Ta, 
which generates particle-hole excitations to all orders, is truncated at the singles and dou¬ 
bles (CCSD) level, Ti -I- T 2 , and includes triples excitations, T 3 , in a nonperturbative but 
approximate way. At the CCSD level, Ti and T 2 are obtained by solving the CC equations, 
which follow from the reference state having no one-particle-one-hole or two-particle-two- 
hole excitations. Equations-of-motion CC methods can access ground and excited states of 
one- or two-particles-attached and one-particle-removed systems from closed-shell nuclei. 

A novel ab initio many-body method is the in-medium similarity renormalization group 
(IM-SRG) II25II26II . The IM-SRG uses a continuous unitary transformation U{s), parame¬ 
terized by the flow parameter s, to drive the Hamiltonian to a band- or block-diagonal 
form. This is accomplished by solving the flow equation = [r;(s), jT(s)], where 

r,(s) = ^C/t(s) is the generator of the transformation. With a suitable choice of ri{s), 
the off-diagonal part of the Hamiltonian is driven to zero as s —^ 00 . Similarly to the CC 
approach, the IM-SRG decouples the closed-shell ground state from the space of particle- 
hole excitations on top of it. Recently, a multi-reference formulation (MR-IM-SRG) was 
developed, which enables to describe also ground states of open-shell nuclei iP?)). 

In self-consistent Green’s function theory (SCGF) (I28II . the quantity of interest is the 
single-particle Green’s function, which describes the propagation of single-particle and 
single-hole excitations in the many-body system. From this the ground-state energy can 
be calculated via the Koltun sum rule II29|I . The Gorkov formalism H.IOI I31II allows to treat 
pairing correlations explicitly and extends SCGF calculations to open-shell nuclei. 

Other calculations of oxygen isotopes in large many-body spaces include the importance- 
truncation no-core shell-model (IT-NCSM) (I32II33II . which extends an exact diagonalization 
with the NCSM by using importance sampling to access larger spaces, and nuclear lattice 
simulations dn, which solve the many-body problem of nucleons on a Euclidean space- 
time lattice. In lattice EFT, the energies of the ground and excited states are obtained by 
propagating the system in imaginary time, as in lattice QCD. 

Approaches based on valence-space Hamiltonians allow the calculation of ground and 
excited states of all nuclei in the valence space provided that the diagonalization is feasi¬ 
ble (|35ll . Many-body perturbation theory (MBPT) (I36II has been used to derive valence- 
space Hamiltonians, based on a diagrammatic approach for calculating the single-particle 
energies (SPEs) and the interactions between valence nucleons. The MBPT includes the 
contributions from configurations outside the valence space perturbatively. State-of-the-art 
MBPT results include NN and 3N forces to third order for the valence-space interactions 
and the consistent SPEs [see Reference (I37II for oxygen isotopes]. 

Recently, nonperturabtive derivations of valence-space Hamiltonians have been achieved 
based on the IM-SRG and CC methods. In the IM-SRG for open-shell nuclei, states with A„ 
particles in the valence space are additionally decoupled from those containing non-valence 
admixtures (1381139II . This gives the energy of the closed-shell core (as in the standard IM- 
SRG), but also valence-space SPEs and interactions. In the CC effective interaction (CCEI) 
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Figure 3 

Top panel: Ground-state energies of oxygen isotopes measured from including experimental 
values of the bound Figure taken from (l42ll . Energies obtained from (a) phenomenological 

forces SDPF-M (|43ll and USDB (|44fl . (b) a G matrix and including Fujita-Miyazawa 3N forces due 
to A excitations, and (c) from low-momentum interactions Viowfc including N^LO 3N forces 
as well as only due to A excitations. The changes due to 3N forces based on A excitations are 
highlighted by the shaded areas. Bottom panels: Left: Ground-state energies of oxygen isotopes 
relative to based on valence-space Hamiltonians, compared to the atomic mass evaluation 
(AME 2012) 1I45II . The MBPT results are performed in an extended - 5 ^/ 7 / 2^372 valence space (I37II 
based on low-momentum NN-I-3N forces, while the IM-SRG dMl and CCEI 1)40 |l results are in the 
sd shell from a SRG-evolved NN-|-3N-full Hamiltonian. Right: Ground-state energies obtained in 
large many-body spaces: MR-IM-SRG l|27|l . IT-NCSM (|27[l. SCGF Il29[l . GC H40|l . based on the 
SRG-evolved NN-|-3N-full Hamiltonian, and Lattice EFT II34II . based on NN-I-3N forces at N^LO. 


approach, a similar decoupling for a given nucleus is achieved using a Lee-Suzuki similarity 
transformation from the CC solution in the large basis space to the valence space dlolllll). 

2.1. Location of the neutron dripline 

The neutron drip line evolves regularly from light to medium-mass nuclei except for a strik¬ 
ing anomaly in the oxygen isotopes, where the dripline is at a doubly magic nucleus and 
anomalously close to the stable nuclei. This anomaly is challenging to explain in microscopic 
theories based only on NN forces that reproduce NN scattering (|42ll. This is illustrated in 
the top panel of Figure]^ with sd-shell calculations based on second-order MBPT [dashed 
lines in panels (b) and (c)], where the ground-state energies decrease up to V = 20, leading 
to an incorrect dripline at ^®0. This is in contrast to phenomenological interactions ad- 
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justed to experiment, shown in panel (a), which have a minimum at because the 0 ( 3/2 
orbital remains unbound. The comparison shows that NN forces only lead to too attractive 
interactions among valence neutrons, which causes the 0 ( 3/2 orbital to become bound in 
and beyond. These deficiencies have been traced to the monopole components of the 
valence-space Hamiltonian, as these parts are amplified with neutron number, and it has 
been argued that 3N forces are important for the monopole components jMli. This is also 
supported by the cutoff dependence of the monopole components with NN forces only 63 - 
First investigations of neutron-rich oxygen isotopes with 3N forces dm) have shown that 
3N forces lead to repulsive interactions between valence neutrons. This is dominated by 
the long-range parts of 3N forces, as highlighted by the results for Fujita-Miyazawa 3N 
forces 63 due to A excitations in the top panel (b) and (c) of Figure This results 
from the interactions of two valence neutrons with any of the nucleons in the core, which 
corresponds to the normal-ordered two-body part of 3N forces, and the repulsive nature 
for valence neutrons is rather general 63 - The dominance of the normal-ordered two- 
body part, which is enhanced by all core nucleons, can be understood based on phase-space 
arguments for normal Fermi systems (Eel) and was verified in explicit calculations II50[ I51II . 
The top panel (c) of Figure]^ shows that N^LO 3N forces, fit to few-nucleon systems only, 
predict the dripline correctly. It is interesting to note that the same 3N forces also lead to 
repulsion in neutron matter, see Section [4.1[ 

This repulsive 3N-force mechanism was confirmed in more recent calculations based 
on large many-body spaces and with improved MBPT and nonperturbative valence-space 
Hamiltonians, as shown in the bottom panels of Figurej^ All results obtained in large many- 
body spaces with a SRG-evolved NN-|-3N-full Hamiltonian (1271129114(111 (see the bottom right 
panel) lead to the correct dripline position at and the different many-body methods 
agree within a few percent. With NN-|-3N-ind forces only (not shown), all oxygen isotopes 
are underbound with respect to experiment and bound up to ^®0. The results obtained from 
ab initio valence-space Hamiltonians are shown in the bottom left panel of Figure The 
nonperturbative calculations are in the sd shell, while the MBPT calculations are performed 
in an extended sdj-j/ 2 'P'i/i valence space. In this case, the range of predictions is broader, but 
the 3N contributions also lead to an increased repulsion with neutron number. The broader 
range here is due to the different Hamiltonians considered (MBPT vs. IM-SRG/CCEI) as 
well as due to different many-body approximations. 


2.2. Spectroscopy 

Next, we discuss excited states in the neutron-rich oxygen isotopes 22,23,24 q^ which have 
been the subject of recent experiments dsn El El El El [53- Figure g compares the 
measured low-lying states with theoretical calculations for each nucleus. We note that the 
spectra of neutron-rich oxygen isotopes have also been calculated in GC theory with phe¬ 
nomenological 3N forces adjusted to the oxygen isotopes (3Nefi) and including the contin¬ 
uum (15811 . with good agreement to experiment. In Figure]^ we show only the CGEI results, 
as they are based on the same SRG-evolved interactions as in the IM-SRG calculations. 

The first oxygen isotope with closed-shell properties for a non-standard magic number, 
has its first 2"*' state at almost twice the energy as those in and In contrast 
to whose closed-shell nature can be qualitatively well described with NN forces due 
to the large separation between the d 3/2 and S 1/2 orbitals, the spectrum of is not well 
reproduced with NN forces: the first 2^ state is below experiment, and the rest of the 
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Figure 4 

Excited-state spectra of 22 , 23 , 24 q imaged on NN-I-3N forces and compared with experiment. Figures 
adapted from (I39II . The MBPT results are performed in an extended sdfj/ 2 Pz /2 space II37II based 
on low-momentum NN-I-3N interactions, while the IM-SRG ll39ll and CCEI g&ll results are in the 
sd shell from the SRG-evolved NN+3N-full Hamiltonian with huj = 20MeV (CCEI and dotted 
IM-SRC) and huj = 24MeV (solid IM-SRC). The dashed lines show the neutron separation energy. 


spectrum is generally too compressed. Only when 3N forces are included, the 2"*" energy 
is in good agreement with experiment. Results with NN+3N forces are shown in Figure]^ 
for MBPT calculations in an extended sdfjf 2 P 3/2 space II37II or with IM-SRG (I39II and 
CCEI (I40II in the sd shell. For the next excited states, the 2+ — 3^ splitting is somewhat 
large in MBPT and results in an inversion of the 3^ and 0^ states and a too-even spacing 
of the other levels. This is not the case for IM-SRG and CCEI, where the level ordering is 
well reproduced and the spacings between states are close to experiment. Remarkably, the 
IM-SRG and CCEI results are within less than 200 keV when the same harmonic-oscillator 
value huj is used. 

Of particular interest is the spectrum of ^^O, which provides a unique test for theory, 
as it simultaneously reflects the features of the doubly magic ^^O and ^'^O. The shell model 
^^O ground state is dominated by one particle in the si /2 orbit, while the two lowest 
excited states are expected to be a single-particle 5/2'^ one-hole excitation, indicative of 
the strength of the ^^O shell closure, and a higher-lying single-particle 3/2^ one-particle 
excitation, reflecting the strength of the ^^O shell closure. The 5/2^ state lies just above 
the neutron decay threshold, while the 3/2^ state resides in the continuum. Reflecting the 
failure to reproduce the shell closure, all calculations with NN forces only predict a too 
low 5/2^ state. With 3N forces included, shown in the middle panel of Figure]^ the 5/2^ 
state is well reproduced in all calculations (15711551 HDli . The position of the 3/2+ state is 
approximately 1.0 MeV too low in MBPT and 1.0 MeV too high with IM-SRG and GCEI, 
which are again in remarkable agreement. The inclusion of the continuum is expected to 
lower the d ^/2 orbital II59II improving the IM-SRG and GCEI results. 

Finally, the right panel of Figure]^ shows the spectrum of with 3N forces in com¬ 
parison with experiment. All calculations find a clear closed-shell signature of a high 2+ 
energy. Moreover, the spacing to the next 1+ state is well reproduced. Note that due to 
the unbound d ^/2 orbital, continuum effects should be included (1581) . 
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Figure 5 

Ground-state energies of and relative to the ground state, and the energy of the first 
excited state in relative to the ground state. The experimental energies are from 
MoNA/NSCL 11601161II and from R3B/LAND II62II . Results are shown for the different many-body 
methods with NN+3N forces as in the bottom panels of Figure]^ MBPT including also residual 
3N forces ll62ll . IM-SRG dSSll, CCEI igOll, MR-IM-SRG SCGF ll29t. and CC lliOt. 


2.3. Beyond the neutron dripline 

Nuclei beyond the neutron dripline play an important role in understanding the behavior 
of extreme neutron-rich systems. The unbound and isotopes are the current ex¬ 
perimental limits in oxygen (|60ll6ni62)l. and the lifetime of makes it the first candidate 
to exhibit two-neutron radioactivity (|63|. A crucial aspect neglected in most of the cal¬ 
culations discussed here is the coupling to the continuum, which has been shown to play 
an important role for the physics of unbound states (|M1) and specifically for the oxygen 
isotopes (ESI). CC SNefif calculations found a typical contribution from continuum coupling 
to be on the order of 200 keV for the unbound states past ^'*0. 

Figure compares the experimental ground-state energies of to theoretical 

predictions with NN-I-3N forces. In addition to the contribution from 3N forces to the 
SPEs and the two-body interactions of valence neutrons, the MBPT results shown EH) also 
include the small contribution from residual three-valence-neutron forces. These become 
more important with increasing neutron number along isotopic chains dSSl, with a repulsive 
contribution from 0.1 — 0.4 MeV for For the ground state of ^®0, Figure shows 

that MBPT, CC, and SCGF agree well with experiment, while the valence-space predictions 
from IM-SRG and CCEI are modestly too bound and unbound, respectively. For the ground 
state of only MBPT gives a result close to experiment. Nevertheless the overbinding 
seen in IM-SRG, and the underbinding obtained in MR-IM-SRG and GGEI (off the scale of 
the plot) are not unreasonable given expected theoretical uncertainties. Better agreement 
is found for the continuum shell model ESI) and in CC 3Neff calculations (|58|). 

As shown in Figure valence-space NN-I-3N calculations consistently predict a low 
first excited 2'^ state in between 1.0 — 1.6 MeV. USDA and USDB interactions give a 
somewhat higher energy of 1.9 MeV and 2.1 MeV, respectively (EH). Experimentally, events 
have been seen at 4 MeV EH. but in calculations the next state lies above 6 MeV due to 
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Figure 6 

Left panel: Ground-state energies of fluorine isotopes measured from compared to AME 
2012 ((45J1. The MBPT calculations H66fl are performed in a proton sd and neutron sdfY/ 2 P 3/2 
extended space based on low-momentum NN-I-3N forces, while the SCGF results 112911 are with the 
SRG-evolved NN+3N-full Hamiltonian. Right panel: Excited-state spectrum of compared 
with experiment II67II and USDB 1144II . The CC results (I41II are obtained in a large many-body 
space based on the optimized chiral NN+3N forces at N^LO. The IM-SRG results II 68 II are 
calculated in the sd shell from the SRG-evolved NN+3N-full Hamiltonian. 


the large d 3/2 — Si /2 gap. Note that the lowest excited states in have been very recently 
measured with high statistics at RIBF/RIKEN. 

2.4. Neighboring open-shell nuclei 

Fluorine isotopes have one more proton, and nuclear forces provide more binding mainly 
due to the tensor part. No predictions based on NN-1-3N forces exist for the dripline in 
fluorine, which lies at least as far as ^^F, but selected isotopes have been calculated from 
SCGF (I29II and MBPT II66II . These are shown in the left panel of Figure]^ Both SCGF and 
MBPT agree well with experiment through ^®F, with MBPT becoming modestly overbound 
beyond. Future experiments are needed to test the predictions beyond ^®F. 

Large-scale CC results for excited states have been presented for 25,26p (pi[7nj| 

with 

3Neff, and for ^2.24p optimized N^LO NN-I-3N interaction II41II . In addition, the 

spectrum of ^^F has been studied recently with the IM-SRG in comparison to new experi¬ 
mental results (|67|). These are shown in the right panel of Figure]^ where we also compare 
to USDB results. Without 3N forces (not shown), the spectrum is much too compressed and 
the ordering of levels is incorrect for both CG and IM-SRG: The first eight excited states lie 
below 2.0 MeV, in clear contrast to experiment. The NN-I-3N results for CG and IM-SRG 
shown in Figureagree well with experiment, but differences are seen due to the different 
starting Hamiltonians. For CC the predicted first 2+ and 1'*’ states are well reproduced, but 
above these the density and ordering of states begins to deviate from tentative experimental 
spin-parity assignments. While the ground-state energy of ^'^F is overbound by 7.7 MeV in 
IM-SRG, the predicted excited-state spectrum is in remarkably good agreement with the 
new experimental measurements, with all excited states below the one-neutron separation 
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threshold less than 200 keV away from experiment. The only exception is the 0'^ state at 
2.6 MeV, also predicted by USDB, to which the experiment is not sensitive. 


3. ISIEUTRON-RICH CALCIUM ISOTOPES 

The calcium isotopes (Z = 20) provide an excellent region to study shell evolution from 
stability towards the neutron dripline. In addition to the standard doubly magic '*‘^Ca and 
"^®Ca, recent pioneering experiments at rare isotope beam facilities have explored new shell 
closures in exotic ®^Ca and ®^Ca dni ESI). The wealth of spectroscopic data on excited 
states and electromagnetic moments and transitions offers excellent tests for nuclear forces. 

The calcium and neighboring isotopic chains lie at the frontier of theoretical calculations 
with 3N forces. Previous studies with phenomenological valence-space interactions iESIItTI 
ESi or beyond-mean-field calculations GSl provide a good description up to ®^Ca, but 
begin to disagree for more exotic isotopes where data was not available. This provides an 
additional motivation for calculations based on chiral EFT interactions. 

The calcium isotopes have been studied with MBPT in an extended pfgQ /2 

valence space based on the same low-momentum NN-I-3N forces as for oxygen (|76l 17711711 
rTHlITgi) . The extended valence space compared to phenomenological interactions in the pf 
shell (13511731174II suggests the need to treat the grg /2 orbital nonperturbatively in MBPT 
calculations based on nuclear forces. Residual 3N forces between valence particles, which 
are important for very neutron-rich nuclei as discussed in Section [2.3| provide very small 
contributions for the calcium isotopes discussed here l|7T | l7 ^ . In addition to valence-space 
calculations, the calcium isotopes have been studied in large many-body spaces, where all 
nucleons are treated as active. The neutron-rich calcium isotopes from (except for 

have been calculated in CC theory with phenomenological 3Neff adjusted to the 
calcium isotopes (I80II81II . In addition, SCGF (I82II and MR-IM-SRG (I83II provide results for 
even calcium isotopes isotopes based on the same SRG-evolved NN-|-3N-full Hamiltonian 
as for the oxygen isotopes. Finally, based on the same Hamiltonian, there are IT-NGSM 
results for the doubly magic ‘^°’'‘*Ga (I5T1) . 

The MBPT and GG results for the ground-state energies are in very good agreement 
with experiment (177117111791180II . while the SGGF (I82II and MR-IM-SRG (I83II both give 
overbound calcium isotopes, suggesting that the SRG-evolved NN-|-3N-full Hamiltonian is 
too attractive for heavier nuclei. Nevertheless this overbinding is systematic for all isotopes, 
and the experimental trend is reasonably reproduced. 

3.1. Shell structure 

The shell evolution along an isotopic chain can be studied with the two-neutron separation 
energy S 2 n, where a significant decrease occurs past a shell closure. Figure [7| compares 
the experimental S 2 n, with calculated MBPT (EZIEI]), GG dHni), SGGF (HH) and MR-IM- 
SRG (1831) results. A first key feature in Figure]^ is the large decrease in S 2 n from N = 28 
to N = 30, a signature of the N = 28 magic number. This shell closure is not reproduced 
with NN forces only (l76l [Mil l3^ . First calculations with 3N forces (TtHIIMIi showed that 3N 
forces are essential for the N = 28 shell closure. This also holds for the improved MBPT 
calculations of References 1(771 EDI, shown in Figure]^ where 3N forces are included to third 
order in MBPT [3N forces were only included to first order in (1751) ]. While the '*®Ga S 2 n 
has not been calculated in GG, the decrease from ‘*®Ga to ®°Ga reproduces experiment very 
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Figure 7 

Two-neutron separation energy S 2 n iri neutron-rich calcium isotopes as a function of neutron 
number N. The experimental energies lESlIZZlED are compared with MBPT predictions GZIED 
based on low-momentum NN-I-3N forces and CC theory with phenomenological SNgfF (I80II . In 
addition, SCGF II82II and MR-IM-SRG II83II results are shown based on a SRG-evolved 
NN+3N-full Hamiltonian. 


well. Also shown in Figure]^ are SCGF II82|I and MR-IM-SRG (I83II results, based on the 
same SRG-evolved NN-|-3N-full Hamiltonian, which somewhat underestimate (SCGF) or 
overestimate (MR-IM-SRG) the decrease in S 2 n past N = 28. 

The flat behavior in S 2 n from = 30 to A’ = 32 was predicted in calculations with 
3N forces dza Ezi), as well as with phenomenological shell-model interactions dill [HI). 
This was recently confirmed with precision Penning-trap mass measurements of ®^’®^Ga at 
TITAN/TRIUMF (l77ll . which found ®^Ca to be 1.74MeV more bound. This represented the 
largest change in the atomic mass evaluation in the last ten years. More recently, the masses 
of ®®’®^Ga were measured in pioneering multi-reflection time-of-flight mass measurement at 
ISOLTRAP/CERN (1711) . The resulting S 2 n values show a decrease from A = 32 to 77 = 34, 
similar to the one past '‘^Ca. This unambiguously establishes the doubly magic character 
of ®^Ca. The shell closure in calcium at A = 32 had already been suggested based on the 
first excited 2^ energy in ®^Ca dUlHSli and nuclear spectroscopy (186118711 . Figure]^ shows 
that the experimental S 2 n from A = 32 to A = 34 are in excellent agreement with the 
MBPT and CC predictions. We also observe that the different calculations with 3N forces 
start to deviate most in ®®Ca. Therefore, a future mass measurement of ®®Ca would provide 
a key test for theory, as well as direct information about the closed-shell nature of ®'*Ca. 
Phenomenological shell-model interactions ifTmiTi]) also provide a good description of S 2 n 
to ®^Ca, while energy-density functionals (jSHl) tend to predict an almost linear trend in S 2 n, 
missing the steeper decrease at A = 28 and A = 32 (E3- 

Another key observable to assess shell evolution is the excitation energy of the first 2+ 
state in the even isotopes. Figure shows the experimental energies compared to the 
available theoretical calculations with NN-I-3N forces: MBPT (17811 and CC (I80p . The agree¬ 
ment of both calculations to experiment is very good, especially in capturing the high 2+ 
energy in ‘*®Ca and ®^Ca, associated with the A = 28 and A = 32 shell closures. Especially 
interesting is the 2"*" energy in the exotic ^'^Ca, recently measured at RIBF/RIKEN (1721) . 
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Figure 8 

Excitation energy of the first 2+ state in the even calcium isotopes as a function of mass 
number A. The MBPT ll79ll and CC results ilsot corresponding to the S 2 n calculations of Figure]^ 
are compared to experiment from 11721 Isa . 


which is also in very good agreement with both MBPT and CC predictions. The relatively 
high energy of this 2"*" state, and the increase compared to other = 34 isotones, suggests 
a shell closure in calcium at Af = 34 lf72ll. 

From Figures and we see that the predictions of shell evolution are consistent 
when considering different observables. In addition, three-point mass differences, defined as 
—[-B(Af -I- 1, -h B(N — 1,Z) — 2B{N, Z)] with (negative) ground-state energy 
B{N, Z), provide another signature for shell evolution, with relative peaks in associated 
with closed shells. When comparing MBPT results dZHl for aI^^ with experiment, the trend 
with neutron number along the calcium chain is well reproduced, with peaks at Af = 28 
and N = 32, and an anomalously low value for ®^Ca, placed between two shell closures 
and dominated by the low-j pi /2 orbital (I90II . 

With present calculations based on NN-)-3N forces, it is difficult to predict the neutron 
dripline because of the very flat behavior of the ground-state energies past ®°Ca found 
in energy-density functional calculations, MBPT, CC, and MR-IM-SRG II8UI1791183[ I91II . 
so that very small interaction contributions can be decisive. In addition, for such weakly 
bound systems, the continuum, currently included only in CC, plays an important role (fSOll. 
and should also be included in all calculations. Moreover, valence-space Hamiltonian for 
more exotic isotopes will have to consider the Ici 5 y 2 and 2si/2 orbitals (|80l I92|l . Improved 
calculations towards the dripline will allow the exploration of a possible A'^ = 40 shell closure 
in ®®Ca, which is the target of ongoing experimental efforts in neighboring isotopes (PI). 
Finally, CC calculations combined with halo EFT have suggested interesting Efimov physics 
around ®®Ca, with the possibility of ®^Ca being a two-neutron halo nucleus dm). 

In addition to the calcium results, SCGF calculations have been performed for the 
ground-state energies of Ar (Z = 18), K (Z = 19), Sc (Z = 21), and Ti (Z = 22) II82|) . 
and the resulting S' 2 „ generally agree well with experiment, apart from an overestimated 
shell gap at A’ = 20 attributed to the NN-|-3N-full Hamiltonian used. In the 

evolution of the splitting between the lowest 1/2"*" and 3/2^ states calculated with SCGF 
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also successfully describes experiment dsn. 

Finally, 3N forces provide similar repulsive contributions to the ground-state energies 
of the proton-rich = 20 isotones, which have been studied in MBPT including the 
isospin-symmetry-breaking parts of chiral EFT interactions dMl). Ground-state properties 
of heavier neutron-rich systems have also been explored in large-scale calculations with 
NN-I-3N forces. First MR-IM-SRG results along the nickel chain {Z = 28) (I83II . based on 
the same NN-|-3N-full Hamiltonian used in calcium, show a similar overbinding, also seen 
in CC calculations for doubly magic nuclei up to ^®^Sn If96ll. 


3.2. Spectroscopy 

The study of shell evolution is complemented by exploring spectroscopic properties of cal¬ 
cium isotopes based on NN-I-3N forces. Spectra beyond 2'^ states based on chiral NN-I-3N 
forces have been obtained with MBPT 1(7^ and with CG 3Nefr for states dominated by 
ground-state particle-hole excitations in (O. Electromagnetic moments and tran¬ 

sitions have also been calculated with MBPT (17611791 IWll . 

The overall agreement of MBPT and CG spectra with experiment is good, with low-lying 
states in odd-mass nuclei consistent with the shell evolution discussed in Section [3.II Some 
deficiencies are present, e.g., with too low 5/2“ states in in MBPT, but generally 

the results are similar to phenomenological interactions ([7^. Low-lying excited states in 
^^Ti, calculated with MBPT II95II . and in neutron-rich 52,54,56rpj^ obtained from CC II80II . 
are also in good agreement with experiment. Future experiments with rare isotopes can test 
MBPT and CC predictions for several excited states in neutron-rich ®^“®®Ca. In addition, 
there are MBPT predictions for the mostly unexplored excited states in proton-rich A'’ = 20 
isotones up to ^®Ni 11951) . 

Electromagnetic moments provide a complementary test of nuclear forces. Figure 
compares the experimental magnetic moments and electric quadrupole moments of neutron- 
rich calcium isotopes with MBPT predictions, including the very recent measurements at 
COLL APS/ISOLDE (1971) . For comparison, the results with KB3G and GXPFIA interac¬ 
tions are also shown. The lower panel of Figure]^ shows the electric quadrupole moments. 
The experimental linear trend from '*^Ca to '*^Ca and ^®Ca to ®^Ca, a signature of the 
filling of the / 7/2 and P 3/2 orbitals, suggests that the quadrupole moments are dominated 
by the single-particle character of the ground states. The MBPT predictions exhibit a very 
good description of the experimental quadrupole moments, in general similar to KB3G and 
GXPFIA interactions and better for ^^Ca. Note that the theoretical results use the same 
neutron effective charges e„ = 0.5e dszj. Electric quadrupole transitions involving excited 
states in ^®“®°Ca have also been calculated in MBPT with the same effective charges dZSl). 
Similar to phenomenology, the agreement with experiment is reasonable, taking into account 
that the measured B{E2) values vary within a factor of 50. 

The upper panel of Figure [^compares the experimental magnetic moments with theo¬ 
retical results obtained with bare g-factors. Calculations for the lighter isotopes '*^’'*®’^®Ca, 
which assume a ^°Ca core, do not reproduce experiment, suggesting the importance of sd- 
shell degrees of freedom for magnetic moments. This is in agreement with the g-factors of 
the 2+ states in (1981199p . as well as calcium isotope shifts (11001) . For '*^Ca MBPT 

reproduces experiment well, and the predictions for neutron-rich '*®Ca and ®^Ca are in very 
good agreement with very recent measurements (^7)). Magnetic moments with phenomeno¬ 
logical interactions are similar to MBPT except for ®^Ca, where the NN-I-3N calculation 
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Figure 9 

Magnetic moments (upper panel) and electric quadrupole moments (lower panel) of the odd 
calcium isotopes as a function of neutron number N. Experiment is compared to MBPT 
calculations with low-momentum NN-I-3N forces in the pfgg /2 valence space II97II . For 
comparison, p/-shell calculations with phenomenological interactions KB3G m and 
GXPFIA Gl are shown. Figure adapted from dig. 


lies between KB3G and GXPFIA. 

The magnetic dipole B(M1) strength in '*®Ga was first calculated with MBPT in Refer¬ 
ence (I26J, with improvements in Reference (l79ll by including 3N forces to third order. The 
experimental peak at 10.2 MeV (fTUTIl is very well reproduced in the MBPT calculation, but 
with too much strength, suggesting a modest quenching of g = 0.9 in the spin g-factor. 
Phenomenological interactions GXPFIA and KB3 also reproduce the concentrated peak 
and strength, but with a more substantial quenching of g = 0.75 (fTOB . 

This quenching is in contrast with the calculations of the magnetic moments in Figure 
where the bare spin g-factors give a good description of experiment. This inconsistency and 
the sensitivities to different interactions in Figure (compared to the smaller spread for 
the quadrupole moments) point to the need for systematic calculations of magnetic moment 
operators in the valence space. Necessary improvements are the inclusion of electromagnetic 
two-body currents (or meson-exchange currents), which are derived in chiral EFT consis¬ 
tently with nuclear forces, as well as controlled calculations of effective operators. Results 
with chiral two-body currents in light nuclei demonstrate that they provide significant con¬ 
tributions to magnetic moments while first applications to medium-mass nuclei have 

focused on Gamow-Teller transitions (14111103ll . 
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4. NEUTROISI-RICH MATTER AND NEUTRON STARS 


The physics of neutron-rich matter covers a wide range of extremes. At very low densities, 
the average interparticle distance is sufficiently large so that details of nuclear forces are 
not resolved and all properties of the system are governed by the large S-wave scattering 
length. In this universal regime, neutron matter shares many properties with cold atomic 
gases close to the unitary limit, which are the subject of active experimental and theoretical 
studies llTOill . At intermediate densities, which are most relevant for nuclei, the properties 
of neutron and nuclear matter are used to guide the development of energy density func¬ 
tionals, in particular to constrain the physics of neutron-rich systems, which are key for 
understanding the synthesis of heavy nuclei in the universe. At higher densities, far beyond 
nuclear densities, the composition and properties of nuclear matter are still unknown. Ex¬ 
otic states of matter containing strange particles or quark matter may be present. On the 
other hand, neutron matter constitutes also a unique laboratory for chiral EFT, because 
all many-body forces are predicted to N^LO, see Section 1.2 This offers the possibility to 


derive reliable constraints based on chiral EFT interactions for the equation of state (EOS) 
of neutron-rich matter in astrophysics, for the symmetry energy and its density dependence, 
and for the structure of neutron stars, but also makes it possible to test the chiral EFT 
power counting and the hierarchy of many-body forces at densities relevant for nuclei. 

The importance of chiral 3N forces for understanding and predicting nuclei has already 
been discussed in Sections and The same 3N forces play an important role for nu¬ 
clear matter. In particular, the saturation of symmetric nuclear matter is driven by 3N 
forces (ITUSinB. While 3N contributions to neutron matter are smaller, they are crucial 
for the EOS of neutron-rich matter, and thus for the symmetry energy and its density 
dependence discussed in Section [4.2[ and for neutron stars in Section [4.3| 


4.1. Neutron matter properties and theoretical uncertainties 

The left panel of Figure [To| shows the energy per particle of neutron matter up to saturation 
density no = 0.16 fm“®. The results are obtained with different many-body methods based 
on chiral EFT interactions, all with the 500 MeV N^LO potential (1201) at the NN level. In 
case where the energy is shown with bands, for MBPT (red lines, cyan and blue band) (11061 
[Ml EMI) and SCGF (EQZI) results, the theoretical uncertainty of the energy is dominated by 
uncertainties in the low-energy couplings ci and ca, which determine the long-range two- 
pion-exchange parts of 3N forces, and not by truncations in the many-body calculation. 
The red lines and blue band show results including contributions from N^LO 3N forces, 
whereas the cyan band includes all 3N and 4N interactions to N®LO. In addition, for the 
blue band the NN potential has been RG-evolved to a low-momentum scale A = 2.0 fm“^. 
We also show CC (|110|l . MBPT of Reference (llll|l . and SCGF (jl07|l results, including 
N^LO 3N forces. These all lie within the overlap of the blue and cyan band (except for 
the lowest density SCGF point, where the zero-temperature extrapolation may be difficult). 
The determination of the d couplings from tvN scattering is consistent with the extraction 
from NN scattering, see, e.g., the discussion in (El, however with large uncertainties. 
Therefore, the d range for the bands in Figure 10 is taken conservatively: at N^LO (red 
lines and blue band), ci = (0.7 — 1.4) GeV“^ and ca = (2.2 — 4.8) GeV“^ (|116l I106|l [with 
a similar d range for SCGF (11071) 1. and at N^LO (cyan band), ci = —(0.75 — 1.13) GeV“^ 
and C 3 = —(4.77 — 5.51) GeV“^ (11171) . 

Figure [To| shows that chiral EFT interactions provide strong constraints on the EOS of 
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Figure 10 

Energy per particle E/N oi neutron matter as a function of density n based on different chiral 
EFT interactions and using different many-body methods The uncertainty bands in the left panel 
show the energy range based on the 500 MeV N^LO NN potential of Reference II20II and including 
N^LO 3N forces in MBPT (red lines) 1I106II or in the SCGF approach (I107II ' as well as including 
all 3N and 4N interactions to N^LO (llOSlfTO^ (cyan band). The blue band shows the results after 
RG-evolution of the NN potential ilTTlIToet . In addition, we show results obtained in CC 
theory iim and in MBPT of Corragio et ah iITTTIi . When bands are given, these are dominated 
by the uncertainties in the ci couplings in 3N forces. Figure adapted from Reference JToell . The 
right panel shows the energy per particle including NN, 3N and 4N forces at N^LO based on 
different N^LO potentials (cyan, magenta, and black bands). The bands include uncertainty 
estimates due to the many-body calculation, the ci couplings, and by varying the 3N/4N cutoffs. 
For details see (11081 fTO^ . For comparison, results are shown at low densities (see also the inset) 
from NLO lattice (I112II and quantum Monte Carlo (QMC) simulations 1I113II . and at nuclear 
densities from variational (APR) I1114II and auxiliary field diffusion Monte Carlo calculations 
(GCR) (I115l> based on 3N potentials adjusted to nuclear matter properties. 


neutron matter, which are consistent among different many-body methods and considering 
variations of the Hamiltonian. The remarkable overlap of the red lines and the blue band 
indicates that neutron matter is, to a good approximation, perturbative for chiral NN 
interactions with A < 500 MeV, see Reference (fTO for details. This has been benchmarked 
by first quantum Monte Carlo calculations with local chiral EFT interactions (nrslfTToli . 
In addition, there are calculations of neutron matter using in-medium chiral perturbation 
theory approaches with similar results (|120l [T2T|l . 

The right panel of Figure 10 shows the first complete N^LO calculation of the neutron 
matter energy, which includes all NN, 3N and 4N interactions to N^LO II108I llUtHl . The 
energy range is based on different NN potentials, a variation of the d couplings (which 
dominates the total uncertainty), a 3N/4N-cutoff variation, and the uncertainty in the 
many-body calculation. We note that the individual 3N topologies at N^LO (see Figure]^ 
give significant contributions to the energy (11081ITO^ . The N^LO range in the right panel 
of Figure]^ is in very good agreement with NLO lattice results ifTT^ and quantum Monte 
Carlo calculations (ESI) at very low densities (see inset), where the properties are deter¬ 
mined by the large scattering length and effective range dSSli. We also find a very good 
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Figure 11 

Energy per particle E/N of neutron matter as a function of density n obtained from quantum 
Monte Carlo calculations with phenomenological NN and 3N potentials i[TT3t . where the strength 
of the short-range 3N interactions is adjusted to the chosen values for the symmetry energy Sy 
(^Sv = 30.5MeV is the symmetry energy obtained for this NN potential.) For each value of Sv, the 
band shows contributions of different 3N operator structures. The inset illustrates a correlation 
between Sv and the L parameter. For details see tmii. 


agreement with other ab initio calculations of neutron matter based on the Argonne NN and 
Urbana 3N potentials: The results based on variational calculations (APR) (11141) are within 
the N^LO band. In addition, the results from auxiliary held diffusion Monte Carlo calcula¬ 
tions (GCR) II115II are shown based on nuclear force models adjusted to an energy difference 
of 32 MeV between neutron matter and the empirical saturation point, see Figure |11[ 
The properties of neutron matter impact the neutron distributions in nuclei. In partic¬ 
ular, a higher neutron matter pressure at typical nuclear densities pushes neutrons further 
out and thus implies larger neutron skins (11231 [HU). Using these correlations the neu¬ 


tron matter results shown in Figure 10 (blue bands in the left panel) predict the neutron 
skin of ^°®Pb to 0.17 ± 0.03fm (fTTO)) . This is in excellent agreement with the extraction 
of 0.156tg ;o 2 ifm from the dipole polarizability dHSI. The theoretical uncertainty is also 
smaller than the target goal of a new PREX measurement using parity violating electron 
scattering at JLAB (11261) . Moreover, including properties of doubly magic nuclei as con¬ 
straints, in addition to low-density neutron matter results, leads to even tighter predictions 
for the neutron skins of ^®®Pb and '*®Ca to be 0.182 ±0.010 fm and 0.173 ± 0.005 fm, respec¬ 
tively dHzl). 


4.2. Symmetry energy 

The symmetry energy characterizes the behavior of the energy of nuclear matter as a func¬ 
tion of proton fraction x = n^/n or asymmetry parameter P = [N — Z^jA = (1 — ‘lx'). 


www.annualreviews.org • Neutron-rich nuclei and neutron-rich matter 19 









Table 1 Predicted range for the symmetry energy Sv and the L parameter, which 
determines the density dependence of the symmetry energy. The results ilTMli are 
obtained from neutron matter calculations based on chiral NN-j-SN forces using the 
expansion Equation!^ with different 7 values, which lead to different incompressibilities 
K. As shown, the predicted ranges for Sv and L depend very weakly on 7 . Also given 
is the predicted range for the L parameter based on quantum Monte Carlo calculations 
with phenomenological NN-[-3N potentials, where the strength of the short-range 3N 
interactions was adjusted to the chosen range for Sv (irm (see also Figure 


7 K [MeV] 

Sv [MeV] 

L [MeV] 

1.2 210 

4/3 236 

1.45 260 

29.7- 32.8 

29.7- 33.2 

30.1 - 33.5 

32.4- 53.4 

32.5 - 57.0 

33.6 - 56.7 

Gandolfi et al. (2011) 

32.0 - 35.1 

40.6 - 63.6 


Around symmetric matter, the energy per particle can be expanded in the following form: 


E{x, n) 

A 


E(0, n) 
A 


+ /9 -E'sym(n) + . . . . 


( 1 ) 


Since the symmetry energy is a key quantity for the equation of state and for astrophysical 
applications, it has been subject of many experimental and theoretical studies 
[T301| . Sometimes the above expansion is truncated after the quadratic term, however, in 
general the symmetry energy Sv is defined as the second derivate with respect to proton 
fraction. 


S, 


1 d'^EIA{n,x) 
8 dx^ 


n = l,x = l/2 


( 2 ) 


where n = n/no is the density in units of the saturation density. Another important 
quantity for astrophysics is the density dependence of the symmetry energy, characterized 
by the L parameter. 


L 


3d^E/A{n,x) 
8 dndx^ 


n = l,x = l/2 


(3) 


To obtain Sv and L it is necessary to extend the microscopic calculations to finite 
proton fractions. This can be achieved by either performing calculations for several different 
asymmetries or by using neutron matter and symmetric nuclear matter as anchor points and 
interpolating these results to arbitrary proton fractions. One such interpolation has been 
used in Reference ilToel) . with an empirical parameterization that includes kinetic energy 
terms plus an interaction energy that is quadratic in the neutron excess 1 — 2x: 


EfA{n,x) 

% 


3 p/3 + (1 _ ^)5/3j (2^)2/3 

— [{2a — 4Qi)a;(l — x) + ol] n + [{2rj — 477n)x(l — a;) + til] fE . 


(4) 


Here, To = 36.84 MeV is the Fermi energy of symmetric nuclear matter at the saturation 
density. The parameters a,rj,aL, and can be determined from the empirical saturation 
properties of symmetric nuclear matter combined with microscopic calculations of neutron 
matter ifToel) . This strategy takes advantage of the fact that the theoretical uncertainties 
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Figure 12 

Symmetry energy as a function of density n obtained from ab initio calculations of asymmetric 
matter based on chiral NN+3N forces H135fl (red band), from microscopic calculations performed 
with a variational approach [Akmal et ah (1998)] 1I114II and from Brueckner-Hartree-Fock 
calculations (BHF) H136II based on the Argonne uis NN and Urbana UIX 3N potentials (with 
parameters adjusted to the empirical saturation point). For comparison, the band over the density 
range n = 0.04 — 0.16 fm~^ is based on an analysis of isobaric analog states (IAS) and including 
the constraints from neutron skins (IAS + skins) (I137II . Figure taken from (11351) . 


of neutron matter calculations are significantly smaller than for systems with finite proton 
fractions. As shown in Table the predicted range for Sv and L depends only weakly 
on the particular choice of 7 , which is correlated with the incompressibility. Also given 
is the predicted range for the L parameter obtained from quantum Monte Carlo calcula¬ 
tions (ITIB based on 3N potentials adjusted to the chosen range of Sv Figure [ll] shows 


these results for neutron matter in more detail, highlighting the correlation between S-u and 
L. It is remarkable how well the Sv — L regions in Table agree, given the very different 
Hamiltonians used in the quantum Monte Carlo calculations. 

The Sv — L predictions based on NN+3N forces agree well with constraints obtained 
from energy density functionals for nuclear masses ifTHTI) and from the ^°®Pb dipole po¬ 
larizability (I125II . In addition, there is good agreement with studies of the Sn neutron 
skin ifTSll . of giant dipole resonances (Ussi), and with an estimate obtained from modeling 
X-ray bursts and quiescent low-mass X-ray binaries A detailed discussion can be 

found in Reference (UMI. 

Recently, the symmetry energy has also been studied in ab initio calculations of asym¬ 
metric matter at small proton fractions based on chiral EFT interactions dSSl). The energy 
of asymmetric matter was found to compare very well with the quadratic expansion even 
for neutron-rich conditions, which was then used to extract the quadratic symmetry energy 
term Esym- In addition, the results were used to benchmark the empirical parameteriza¬ 
tion, Equation]^ with very good agreement. In contrast to other calculations, the results 
are based on 3N forces fit only to light nuclei, without adjustments to empirical nuclear 
matter properties. The results for Esym are compared in Figure |l2| with constraints from 
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a recent analysis of isobaric analog states (IAS) and including the constraints from neu¬ 
tron skins (IAS -I- skins) (I137D . showing a remarkable agreement over the entire density 
range. Compared to extracting the symmetry energy from neutron matter calculations 
using the empirical parameterization, Equation the uncertainty is reduced due to the 
explicit information from asymmetric matter. 


4.3. Neutron stars 


The structure of nonrotating neutron stars can be studied by solving the Tolman- 
Oppenheimer-Volkov equations (11381) . based on an EOS of neutron star matter, i.e., matter 
in beta-equilibrium. The proton fraction x for matter under these conditions is determined 
by minimizing, for a given nucleon density, the total energy per particle plus the contribu¬ 
tions from electrons and from the rest mass of the nucleons. The resulting proton fractions 
are of the order x « 5% (iTTHinro . 

Since the central densities of neutron stars can significantly exceed the regime for which 
reliable microscopic nuclear matter calculations are possible, it is necessary to extend the 
EOS systematically to higher densities. This can be achieved by using microscopic results up 
to a maximal density and employing general piecewise polytropic extensions beyond (HH 
ITMl) . This strategy allows for soft regions and generates a complete set of possible EOSs 
at high densities, independent of the assumptions on the interactions and constituents of 
matter at high densities. In the end only those EOS are retained that (a) remain causal for 
all relevant densities, and (b) are able to support a neutron star mass M = M, the mass of 
the heaviest observed neutron star. 

The left panel of Figure |13| shows the resulting uncertainty band for the pressure 
as a function of mass density. The blue band at lower densities represents the pressure 
predicted for matter in beta equilibrium based on chiral EFT interactions (see red dashed 


lines in the left panel of Figure 10). The bands at higher densities give the EOS range. 


which is the envelope of all allowed polytropes at higher densities. The lighter blue band 
at high densities corresponds to the mass constraint M = I.QTMq, the central value of 
the two-solar-mass neutron star measured by Shapiro delay and the lower Icr mass 

of the recently observed most massive neutron star from radio timing observations (fliOl) . 
whereas the darker blue band corresponds to M = 2.4Mq, a fictitious heavier neutron 
star. Obviously, the higher the mass of the heaviest neutron star observed, the stronger the 
EOS band is constrained. This uncertainty band is compared with a representative set of 
EOSs used in the literature. This set contains EOSs calculated within diverse theoretical 
approaches and based on different degrees of freedom. For details and notation we refer 
to Reference ifml) . We find that a significant number of EOSs are not compatible with 
the lower density band based on chiral EFT interactions. In addition, at higher densities 
only very few EOSs, including the variational EOSs based on phenomenological nuclear 
potentials dm AP3 and AP4 in the left panel of Figure |13[ are within the uncertainty 
bands over the entire density range. Finally, these constraints imply that a 1.4 (1.97) Mq 
neutron star does not exceed densities beyond 4.4 (7.6) no, which corresponds to a Fermi 
momentum of only 550 (660) MeV. 

The EOS bands directly translate into constraints for the radii of neutron stars. In the 
right panel of Figure we present the radius constraints obtained from the EOS bands 
shown in the left panel for the case M = 1.97 Mq (blue dashed lines) and the corresponding 


bands based on complete N'^LO calculations shown in the right panel of Figure 10 (red 
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Figure 13 

Left panel: Constraints for the pressure P of neutron star matter as a function of mass density p 
compared to EOSs commonly used to model neutron stars II141II . The blue band at low densities 
represents the pressure predicted by the neutron matter results of Figure [To| (red lines of the left 
panel) and incorporating beta equilibrium. The bands at high densities are the envelope of general 
polytropic extensions that are causal and support a neutron star of mass M = 1.97 Mq (light blue 
band) and M = 2.4Mq (darker blue band). For details see II1Q6D . Right panel: Constraint s on the 
mass-radius diagram of neutron stars based on the N^LO neutron matter results of Figure [To| 
(right panel) and following |[T06ll for the extension to neutron-star matter and to high densities 
(red band), in comparison to the constraints based on RG-evolved NN interactions and 3N forces 
at N^LO (thick dashed blue lines, based on the 1.97 Mq band of the left panel). We also show the 
mass-radius relations obtained from EOSs for supernova simulations. For details see lfT09ll . 


band). The radius uncertainty band represents an envelope of a large number of individual 
EOS reflecting the uncertainties at low densities and in the polytropic extensions to high 
densities 111061 llOQll .The results of Figure 13 predict a radius range of 9.7 — 13.9km for 
a IAMq neutron star based on the complete N^LO calculations. The largest supported 
neutron-star mass is found to be 3.1 Mq, with a corresponding radius of about 14km. 
These results agree very well with the mass-radius constraints from the neutron matter 
calculations based on RG-evolved NN interactions with N^LO 3N forces. Furthermore, the 
radius constraints are also consistent with astrophysical extractions obtained from modeling 
X-ray burst sources, see, e.g.. Reference ifmli . 

For astrophysical applications it is crucial to reduce the theoretical uncertainties of the 
EOS over the entire density range shown in the left panel of Figure [T^ At lower densities this 
requires improved estimates and a reduction of the theoretical uncertainties in nuclear forces. 
Specifically, this involves the inclusion of higher-order contributions in the chiral expansion, 
systematic order-by-order convergence studies, as well as improved determinations of the 
low-energy couplings. At higher densities, novel observations are expected to provide new 
and model-independent constraints: On the one hand, the observation of heavier neutron 
stars leads to a systematic reduction of the uncertainties, as illustrated in the left panel of 
Figure On the other hand, astrophysical information on neutron star radii will provide 
significant constraints on the EOS. In particular, the gravitational wave signal from mergers 
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of binary neutron stars II142111431114411145t and neutron star-black hole mergers (I146II has 
been shown to be sensitive to properties of the EOS at high densities and to neutron star 
radii. Hence, a detection with, e.g., advanced LIGO (1T47I) will significantly improve our 
knowledge of the EOS in this regime. Moreover, future observations with the Neutron star 
Interior Composition ExploreR (NICER) and the Large Observatory for X-Ray Timing 
(LOFT) offer completely new perspectives for measuring neutron star radii. This will in 
turn provide important insights into nuclear forces at neutron-rich extremes. 


5. SUMMARY AND OUTLOOK 

We have shown that the physics of 3N forces connects neutron-rich nuclei with neutron-rich 
matter in neutron stars. The main features discussed in this review can be summarized by 
3N forces having two effects for neutron-rich systems: first, 3N forces provide a repulsive 
central interaction, which drives saturation and leads to a stiffening of the neutron matter 
equation of state with increasing density. The same repulsion is important for the contribu¬ 
tions from 3N forces to the ground-state energies that are key for the location of the neutron 
dripline. Second, 3N forces provide an attractive spin-orbit interaction, which increases the 
spin-orbit splittings and gives rise to the associated magic numbers, e.g., at N = 28. This 
physics and its interplay with the repulsive central interaction is at play for the formation 
and evolution of shell structure and for the spectroscopy of neutron-rich nuclei discussed 
in this review. The experimental discoveries for neutron-rich oxygen and calcium isotopes 
show that 3N forces provide an exciting link between the theoretical frontier in effective field 
theories and many-body methods with the exploration of exotic nuclei at rare isotope beam 
facilities worldwide. Three-nucleon forces are also key for the properties of neutron-rich 
matter at nuclear densities, which impacts the neutron skin, the symmetry energy, and the 
structure of neutron stars, in particular their radii. 

Finally, we list a number of open problems and opportunities. First, it is important 
to study the order-by-order convergence in chiral EFT, and to carry out first complete 
N^LO calculations of nuclei. This is especially important, because the contributions are 
known to be important for an accurate reproduction of NN phase shifts at energies relevant 
to nuclei. It is crucial to study the theoretical uncertainties due to the truncation in chiral 
EFT, due to uncertainties in the low-energy couplings in NN and 3N forces, and due to the 
many-body calculation. While the latter has been well documented for many approaches 
and with benchmarks (see, e.g., the remarkable agreement for the same Hamiltonian with 
different ab initio methods discussed here), the theoretical uncertainties in the Hamiltonian 
have been less explored. Because most calculations so far are based on the N^LO potential 
of Reference l|20ll. it is also necessary to explore different NN potentials and many new 
developments in chiral EFT, including the optimized potentials of References dm fliSll . 
new local chiral potentials (nm, and improved chiral potentials of Reference 1 (1491) . all 
with corresponding 3N forces. In addition, it is important to explore improved power 
counting schemes (e.g., see the discussion in Reference (UMI) and to investigate nuclear 
forces with explicit A degrees of freedom (11511 fTS^ . Moreover, the discussed results show 
that additional work is needed to quantify the theoretical uncertainties, especially due to 
the truncation in chiral EFT (IT^ . which is particularly relevant to the long-range parts 
of 3N forces due to the large Ci values. 

Important advances for ab initio calculations are the investigation of open-shell nuclei 
with a range of many-body methods, the inclusion of the continuum for loosely bound 
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and resonant states dose to the dripline (El [SHI, the study of electroweak transitions with 
effective operators and consistent two-body currents based on chiral EFT (I154II10J1|1U2[HT1| . 
extending the calculations to heavier nuclei with a controlled convergence in terms of 3N 
matrix elements included, and the reduction of uncertainties in the equation of state to 
further constrain the properties of neutron-rich matter and neutron-stars. Finally, it would 
be very interesting to transport the knowledge and constraints from nuclear forces to density 
functional calculations of all nuclei, especially regarding the structures of 3N forces and their 
impact on neutron-rich systems. 
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